Entropy dependence of correlations in one-dimensional SU(N) antiferromagnets 



Laura Messio 1 and Frederic Mila 2 

'institut de Physique Theorique (IPhT), CEA, CNRS, URA 2306, F-91191 Gif-sur-Yvette, France 
2 Institute of Theoretical Physics, Ecole Poly technique Federale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland 

(Dated: July 6, 2012) 

Motivated by the possibility to load multi-color fermionic atoms in optical lattices, we study the entropy 
dependence of the properties of the one-dimensional antiferromagnetic SU(N) Heisenberg model, the effective 
model of the SU(N) Hubbard model with one particle per site (filling 1/N). Using continuous-time world 
line Monte Carlo simulations for N = 2 to 5, we show that characteristic short-range correlations develop at 
low temperature as a precursor of the ground state algebraic correlations. We also calculate the entropy as a 
function of temperature, and we show that the first sign of short-range order appears at an entropy per particle 
that increases with N and already reaches 0.8ks at N = 4, in the range of experimentally accessible values. 
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Lattice SU(N) models play an ever increasing role in 
the investigation of strongly correlated systems, both in con- 
densed matter and in cold atoms. The first systematic use of 
these models took place in the context of the large- N gen- 
eralization of the SU(2) Heisenberg model, in which con- 
jugate (or self-conjugate) representations are put on the two 
sublattices of the square lattice so that a SU(N) singlet can 
be formed on two sites " . Over the years, another class of 
SU(N) models with the same representation at each site has 
appeared as the relevant description of the low temperature 
properties in several contexts. In particular, the SU(S) model 
corresponds to the spin-1 Heisenberg model with equal bilin- 
ear and biquadratic interactions ~ , while the SU (4) model is 
equivalent to the symmetric version of the Kugel-Khomskii 
model of Mott insulators with orbital degeneracy ' . These 
models have however attracted renewed attention recently 
as the appropriate low energy theory of ultracold gases of 
alkaline-earth-metal atoms in optical lattices in the Mott in- 
sulating phase with one atom per site, the parameter N corre- 
sponding to the number of internal degrees of freedom of the 
atoms 9 . 

A peculiar characteristic of these SU (N) models is that one 
needs N sites to form a singlet. This is often reflected in their 
ground state properties. In one dimension, the SU (N) model 
has been solved with Bethe ansatz for arbitrary N , and the 
dispersion of the elementary fractional excitations has a peri- 
odicity 2tt/N. On a ladder, the SU(4) model has a plaque- 
tte ground state 11 . In two-dimensions, the SU(3) model on 
both the square and triangular lattices has long-range color 
order with 3 -site periodicity along the lines ' , while on the 
kagome lattice it is spontaneaously trimerized . The SU (4) 
model on the checkerboard lattice also has a plaquette ground 
state 2 . Even on the square lattice, where the ST/" (4) model 
undergoes spontaneous dimerization with possibly algebraic 
correlations , neighboring dimers involved pairs of differ- 
ent colors, so that the 4 colors are indeed present with equal 
weight on all plaquettes. The general properties for arbitrary 
N are not known however. An adaptation of the previous 
large- N studies has been proposed for m atoms per site 15 . If 
m = O(N), the ground state has been proposed to be a chiral 
spin liquid for large N. 

The wealth of ground states predicted for different N on 



various lattices calls for an experimental investigation. Ultra- 
cold fermionic atoms can a priori lead to very accurate realiza- 
tions of these models. However, the temperature is a limiting 
factor. It can be lowered with respect to the initial temper- 
ature if the optical lattice is adiabatically switched on 16 , but 
it cannot be made arbitrarily small. In fact, with adiabatic 
switching, one can control the entropy rather than the tem- 
perature, and in current state-of-the-art experimental setups, 
the lower limit for fermions with N — 2 is equal to 0. 77 ks 
per particle 17 . If contact is to be made with experiments on 
cold atoms, it is thus crucial to know the properties of a given 
model as a function of entropy. For the SU(2) Heisenberg 
model on the cubic lattice, Neel ordering takes place at an en- 
tropy 0.338 ks, i.e. about half the value that can be achieved 
today . 

The first hint that increasing the number of colors might 
help in beating this experimental limit has been obtained in the 
context of a high temperature investigation of the Af-flavour 
Hubbard model by Hazzard et al , who have shown that the 
effective temperature reached after introducing the optical lat- 
tice decreases with N under fairly general conditions. How- 
ever, to the best of our knowledge, no attempt has been made 
so far to determine how the temperature or the entropy below 
which signatures of the ordering will show up depends on N. 

In this Letter, we address this issue in the context of the 
one-dimensional (ID) antiferromagnetic SU(N) Heisenberg 
model on the basis of extensive Quantum Monte Carlo (QMC) 
simulations. As we shall see, the ground state algebraic corre- 
lations lead to characteristic anomalies in the structure factor 
upon lowering the temperature. These anomalies only become 
visible at quite low temperature, but remarkably enough, the 
corresponding entropy per particle increases with N, leading 
to observable qualitative effects with current experimental se- 
tups for N >4. 

The SU (N) Heisenberg model — A good starting point to 
discuss -color fermionic atoms loaded in an optical lattice 
is the SU (N) Hubbard model defined by the Hamiltonian: 

H = t ^ $ ai Caj + h.c.) + U ^2 noting, (1) 

(ijj)a i,a<[3 

where cf a and q jQ! are creation and annihilation operators of 
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a fermion of color a = 1 . . . N on site z and the sum is over 
the first-neighbors of a periodic chain of length L. n a i is the 
number of fermions of color a on site i. At filling 1/N, i.e. 
with one fermion per site, the ground state is a Mott insulator, 
and to second order in t/U, the low-energy effective Hamil- 
tonian is the SU (N) Heisenberg model with the fundamental 
SU(N) representation at each site, and with coupling con- 
stant J = 2t 2 /U. Setting the energy unit by J = 1, this 
Hamiltonian can be written (up to an additive constant): 

H = J2 p ir (2) 

where permutes the colors on sites i and j. If we denote 
by the operator that replaces color ft by a on site i, this 
permutation operator can be written as: 

This effective Hamiltonian is an accurate description of the 
system provided the temperature is much smaller than the 
Mott gap. In terms of entropy, the criterion is actually quite 
simple. The high temperature limit of the entropy per site of 
the SU (N) Hubbard model at 1 /TV-filling can be shown to be 
equal to k B (N In N - (N - 1) ln(N - 1)), while that of the 
SU(N) Heisenberg model is equal to fc# In N. So we expect 
the description in terms of the Heisenberg model to be accu- 
rate when the entropy is below kslnN. For SU (2), this is a 
severe restriction for experiments since In 2 ~ 0.693..., but al- 
ready for SU (3), this is less of a problem since In 3 ^ 1.099. 
Of course, this is not the whole story since what really matters 
is the entropy below which specific correlations develop, but 
this is an additional motivation to consider SU (N) models 
with TV > 2. 

Exact results. — A number of exact results that have been 
obtained over the years on the ID SU (N) Heisenberg model 
will prove to be useful. The model has been solved with Bethe 
ansatz by Sutherland . He showed that, in the thermody- 
namic limit, the energy per site is given by 

k=2 

where £ is the Riemann's zeta function. Some values are given 
in Tab. II. In addition, he showed that there are N — 1 branches 
of elementary excitations which all have the same velocity 
v = 2tt/N at small k. Affleck has argued that the central 
charge c should be equal to N — 1 , and Lee has shown that 20 , 
at low temperature T, the entropy is given by: 

S(T) = kBN{N - 1) T + 0(T% (5) 
6 

a direct consequence of c = N — 1 and v = 2n/N since the 
linear coefficient is equal to ttc/3v. 

The QMC algorithm. — Quantum Monte-Carlo is the most 
efficient method to study the finite temperature properties of 
interacting systems provided one can find a basis where there 



is no minus sign problem, i.e. a basis in which all off-diagonal 
matrix elements of the Hamiltonian are non-positive. For the 
SU (2) antiferromagnetic Heisenberg model on bipartite lat- 
tices, this is easily achieved by a spin-rotation by tt on one 
sublattice. For SU(N) with N > 2, there is no such gen- 
eral solution, but in ID one can get rid of the minus sign on 
a chain with open boundary conditions, as already noticed for 
the SU(4) model . Let us start from the natural basis con- 
sisting of the N L product states 0i|a^) = |ao, . . . , 
where is the color at site i. In this basis, all off-diagonal 
elements of the SU (N) model of Eq.2 are either zero or pos- 
itive. However, a generalization of the Jordan- Wigner trans- 
formation allows to change all these signs on an open chain. 
This transformation is defined by: 

|ao, - - - , (-l)K«o,..,^-i) | ao? 5 aL _ l)? (6) 

where r(ao? • • • ? a L-i) is the number of permutations be- 
tween different color particules on neighboring sites needed to 
obtain a state such that the are ordered (c^ < ctj fori < j). 
This basis change is equivalent to a Hamiltonian transforma- 
tion, the new Hamiltonian being given by: 

H = Y,Y1 \ SrSj a -^Sfsf a \ . (7) 

(ij) a \ [3^a J 

On a periodic chain, the equivalence with the Hamiltonian of 
Eq. 2 is not exact, but the difference disappears in the ther- 
modynamic limit. So in the following we will simulate the 
Hamiltonian of Eq.7. 

To do so, we have developed a continuous time world-line 
algorithm with cluster updates adapted to the model of Eq. 7 
with N colors. The partition function Z is expressed as a path 
integral over the configurations (j) : r —> <p(r), where r is 
the imaginary time going from to f3 = and <p(r) is a 
basis state. The functions <j) that contribute to the integral can 
be represented by 0(0) and by a set of world line crossings 
{(i, j, t)} that exchange the colors of two sites i and j at time 
r. A local configuration c on a link ij at time r is represented 
by 

fa i {T+)a j (T+)\ 
C ~ { ai (r-) aj (r-)) (8) 

Cluster algorithms are well documented for 2-color models. 
Here we generalize the approach to SU (N) by choosing ran- 
domly two different colors p and q out of and by construct- 
ing clusters on which only these two colors are encountered. 
The steps to construct the clusters are the following. We first 
randomly place elementary graphs in the configuration using 
a Poisson distribution. These graphs are drawn in the first col- 
umn of Tab.I and the Poisson time constant is given in the last 
column. They are accepted only if Aq(c) = 1 (if a color 
which is neither p nor q appears in the local configuration, the 
graph is rejected). Then we assign graphs to the world-line 
crossings between p and q colors using the last two columns 
of Tab.I. At the places where no graph has been attributed, 
we follow the path with the same color. Finally we follow 
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TABLE I: Fabrication rules for clusters. In the first column are 
drawn all possible graphs G. The three following columns are Ac (c) 
equal to 1 when G is compatible with the local configuration c. The 
last column is the constant of time used to sample the graphs with a 
Poisson law. e is a small number ensuring the ergodicity (e is taken 
0.01 in our simulations). 




FIG. 1 : Evolution of the energy per site E and of the entropy per 
site S as a function of the temperature T for different N on a L — 60 
chain. The inset shows the slope of the entropy at T = 0, given in 
Eq. 5. The curvature being positive at T = 0, the curves go higher 
than the tangent (dashed line). 

each constructed cluster and exchange p and q on it with a 
probablity 1/2 (Swendsen-Wang algorithm). This constitutes 
a Monte Carlo step. 

Using this algorithm, we have calculated the energy per site 
E, which is given by: 




where n is the number of Monte Carlo steps and n{(j)) the 
number of world-line crossings in the configuration cj>, the di- 
agonal correlations defined by 



CU) = {Y,SS a Sf a j - ^ (10) 
and the associated structure factor defined by 

^^^lE^'' (id 

3 

This structure factor is normalized in such a way that 

The results. — We have studied chains of length L = 60 
for T from 0.01 to 20 with a number of colors N = 2,3,4 
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BA(L = 60) 
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-0.8251193 




-0.82577(2) 
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-0.884730 




-0.88541(2) 



TABLE II: Ground state energies per site obtained for several N 
with the Bethe Ansatz (BA) in the thermodynamic limit and on a 
finite size chain , and at T — 0.01 with the quantum Monte-Carlo 
algortihm (see text) on a L — 60 chain and n — 10 7 Monte-Carlo 
steps. 



and 5, and a number of Monte Carlo steps n at least equal to 
10 6 . The correlation time measured by the binning method 
indicates that around N steps are needed to obtain uncorre- 
cted configurations, whatever the temperature, and that the 
precision on the energy per site E is better than 10 -4 . This 
could be confirmed by the comparison of the limit of the en- 
ergy when T — >> with the exact finite L value for SU (3). 23 
Moreover, the energy of the ground state differs from that of 
the thermodynamic limit by less than 8.10 -4 . So, for our pur- 
pose, the finite size effects can be considered to be negligible 
(see Tab. II). The entropy per site S has been deduced from 
the energy E by an integration from high temperature: 

S(T) = 5(oo)- [°° dr 1 ^^ (12) 
Jt t dr 

where S(oo) = fcsln(iV). E and S are plotted in Fig. 1 
for different N as a function of T. Since the entropy is the 
result of a numerical integration, it is important to check its 
accuracy, especially at low temperature since by construction 
it has to be correct at high temperature. Now, we know that, at 
low temperature, the entropy must be linear with a slope equal 
to ksN(N — l)/6 (see Eq.5). This is confirmed by the inset 
of Fig. lb), in which one clearly sees that the entropies times 
6/ksN(N — 1) lie on top of eachother at low temperature. 

Now, the stabilization of the energy at low T occurs at 
a temperature that decreases when N increases. Thus, one 
could naively think that it will be more difficult to observe 
the development of the ground state correlations when TV in- 
creases. However, this is not true if one considers the en- 
tropy. Indeed, the entropy grows much faster at low temper- 
ature when increasing N. So, the temperature corresponding 
to a given entropy decreases very fast when N increases. 

We now look at the diagonal correlations C(k). They have 
been calculated for different temperatures, but, in view of the 
implications for ultracold fermionic gas, we represent them as 
a function of the entropy per site S. Since the system is ID, 
there is no long range order, hence no Bragg peaks. Never- 
theless, short-range correlation develop at low entropy. They 
translate into finite height peaks in C(k) at finite temperature, 
and singularities at zero temperature. The number and the 
position of these peaks depend on the number of colors N. 
From the Bethe ansatz solution, singularities are expected to 
occur at k = 2pir/N with p = 1, N — 1. The results of 
Fig. 2 agree with this prediction: there is a single peak at 7r for 
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FIG. 2: Evolution of the structure factor C(k) as a function of the entropy per site S for different TV on a L = 60 chain. 
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FIG. 3: Structures factor at low temperature (fcgT = 0.01) for 
different TV on a L = 60 chain, with n = 10 7 Monte-Carlo steps. 
Small peaks are clearly visible at k = tt for TV = 4 and = 4tt/5 
and 67r/5 for N = 5. The data for TV = 4 are in perfect agreement 
with those of Ref .2 1 . 



SU (2), while TV — 1 peaks are indeed present for SU (TV) at 
sufficiently small entropy. Note however that all peaks do not 
have the same amplitude for TV > 4. For TV = 4 and 5, two 
types of peaks not related by the symmetry k —> 2tt — k are 
present. The peaks at 2tt/N and 2 (TV — 1)tt/N are much more 
prominent, and they start to be visible at much larger entropy. 

At the maximal entropy, the structure factor C(k) is flat (see 
Fig. 2). At large but finite entropy, it presents a broad max- 
imum at k = tt for all TV. This reflects the simple fact that 
colors tend to be different on neighboring sites. More specific 
correlations appear upon lowering the entropy. For SU(2), the 
peak at k = tt just gets more pronounced. To observe the de- 
velopment of the singularity typical of the SU{2) ground state 
algebraic correlations will however require to reach rather low 
entropy. This should be contrasted with the TV > 2 cases, 
where a qualitative change in the structure factor occurs upon 



reducing the entropy: the broad peak at k = tt is replaced 
by peaks at 2tt/N and 2(N — 1)tt/N. One can in principle 
read off the corresponding entropy from Fig. 2. To come up 
with a quantitative estimate, we note that, upon reducing the 
entropy, the curvature of the structure factor at k = tt changes 
sign from positive at high temperature to negative when the 
peaks at 27r/TV and 2 (TV — 1)tt/N appear . This occurs at 
S c /k B = 0.58, 0.87 and 1.08 for N = 3, 4 and 5 respec- 
tively. This characteristic entropy S c increases more or less 
linearly with N as S c ~ 0.2Nk B , and for TV = 4 and 5, it lies 
in the experimentally accessible range. This is mostly a con- 
sequence of the temperature dependence of the entropy, which 
grows much faster with TV at low temperature. The character- 
istic temperature at which deviations from the broad peak at 
k = tt occur depends only weakly on TV. Finally, secondary 
peaks appear at lower temperature (see Fig. 3) 

Conclusions. — We have shown that the entropy at which 
the periodicity characteristic of the zero temperature algebraic 
order of SU (TV) chains is revealed increases significantly with 
TV. For TV = 4, this entropy is already larger than the entropy 
per particle recently achieved in the TV = 2 case in the center 
of the Mott insulating cloud (0.77 k B ) • Whether a similar 
entropy can be achieved for TV > 3 remains to be seen. As 
shown by Hazzard et al , if the initial temperature is fixed, 
the initial entropy in a 3D trap increases with TV as TV 1 / 3 , im- 
plying that one might have to go to values of TV larger than 4 to 
reach a final entropy low enough to observe characteristic cor- 
relations. However, evaporative cooling might allow to reach 
initial entropies that are less dependent on TV. In a recent ex- 
periment on 173 Yb, the initial entropy reported by Sugawa et 
al 24 for this TV = 6 case is not much higher than in TV = 2 
experiments . It is our hope that the present results will en- 
courage the experimental investigation of the 1 /TV-filled Mott 
phase of TV-color ultracold fermionic atoms. 
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